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Abstract 

The equations of motion of compact binary systems and their associated Lagrangian formulation 
have been derived in previous works at the third post-Newtonian (3PN) approximation of general 
relativity in harmonic coordinates. In the present work we investigate the binary's relative dynam- 
ics in the center-of-mass frame (center of mass located at the origin of the coordinates). We obtain 
the 3PN-accurate expressions of the center-of-mass positions and equations of the relative binary 
motion. We show that the equations derive from a Lagrangian (neglecting the radiation reaction), 
from which we deduce the conserved center-of-mass energy and angular momentum at the 3PN 
order. The harmonic-coordinates center-of-mass Lagrangian is equivalent, via a contact transfor- 
mation of the particles' variables, to the center-of-mass Hamiltonian in ADM coordinates that is 
known from the post-Newtonian ADM-Hamiltonian formalism. As an application we investigate 
the dynamical stability of circular binary orbits at the 3PN order. 
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I. INTRODUCTION 



The problem of the dynamics of two compact bodies is part of a larger program aimed at 
unraveling the information contained in the gravitational-wave signals emitted by inspiralling 
and/or coalescing compact binaries (see Refs. [l|,li for reviews). The current treatment of 
the problem is Post-Newtonian (expansion when the speed of light c +00; following the 
standard practice we say that a term of order relative to the Newtonian force belongs 
to the nPN approximation). The first breakthrough in the problem of dynamics has been 
;he completion of the equations of motion of two point-like particles up to the 2PN order 

n recent years the (quite involved ) eq uations of motion at the next 
3PN order have also been successfully derived [H Q, EIH El El H El H Q • 

Up to the order 3.5PN there is a clean separation between the conservative part of the 
dynamics, made of the "even" Newtonian, IPN, 2PN and 3PN approximations — with the 
2PN and especially the 3PN ones being very difficult to obtain — , and the part associated 
with the radiation reaction, and consisting of the "odd" 2.5PN and 3.5PN orders (which are 
comparatively much simpler to control than 2PN and 3PN). In principle the conservative part 
of the equations yields a point of dynamical general-relativistic instability at which there is 
(presumably) a transition from the adiabatic inspiral to the final plunge and coalescence. On 
the other hand, the non-conservative terms — i.e. 2.5PN computed in Refs. HQ, 0,0,0,01, 
3.5PN and 4.5PN computed in Refs. I21L I22I |23| — are determined by the boundary 



conditions imposed on the gravitational field at infinity. 

One should not confuse the latter nomenclature for post-Newtonian orders with a different 
one applied to the gravitational field at future null infinity. There the "Newtonian" order, 
which has a quadrupolar wave pattern, corresponds to the dominant odd term in the local 
equations of motion, i.e. 2.5PN, while the IPN order, which is both quadrupolar and 
octupolar, corresponds to 3.5PN in the local equations. And so on. Because of the presence 
of tails at the 1.5PN order in the wave field at infinity there is a contribution at the 4PN 
order in the equations of motion that is "odd" in the sense of being associated with radiation- 
reaction effects 0. Similarly one expects that the known ta,ls-of-ta,ls arising at the 
3PN order in the wave field will correspond to an odd contribution at the 5.5PN order in 
the equations of motion^. 

Two different methods, relying on two independent frameworks, have been applied to the 
equations of motion at the 3PN order. Jaranowski and Schafer Il2j | , working within the 
ADM-Hamiltonian formalism of general relativity, derived the Hamiltonian describing the 
motion of two compact bodies, in ADM coordinates, and in the center-of-mass frame. The 



This difference by 2.5PN orders explains why the equations of motion are insufhcient as regards the 
radiative aspects of the problem. For analyzing the waves emitted by inspiralling compact binaries one 
needs not only the solution of the problem of motion but also the (equally crucial) solution of the problem 
of gravitational- wave generation [3| ■ 
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iamiltonian was later generalized to an arbitrary frame in Ref. Blanchet and Faye 

(following the method proposed in Ref. [^) performed a direct iteration of 
the equations of motion in harmonic coordinates, and in a general frame. The Lagrangian 
of the motion was then deduced from the equations of motion The end results provided 
by these two methods — ADM-Hamiltonian and harmonic-coordinates — have been shown 
to be physically equ.valeut Q Q : there exists a unique -coutaef trausformatiou of 



the binary's dynamical variables that transforms the harmonic-coordinates Lagrangian |2( 
into a different Lagrangian, whose Legendre transform agrees with the ADM-coordinates 
Hamiltonian 



In the works 
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UJ, DJ, UJ and |ia mpact objects are modelled by 

structureless (non-spinning) point-particles. Such a modelling is quite efficient and physically 
sound when describing the inspiral of compact binaries, but the shortcoming is the necessity 



of a regularization for removing the infinite self field of each of the point-particles, 
regularization of Hadamard (or, more precisely, a refined form of it proposed in Refs. 



'he 



19| and implemented in the harmonic-coordinates approach) has been applied but turned 



out to be incomplete in the sense that one (and only one) numerical coefficient remains 
undetermined at the 3PN order : a;, in the ADM-Hamiltonian formalism A in 

he harmonic-coordinates approach jlfil . 0|. This coefficient has been computed in Ref. 

15! with the help of a dimensional regularization instead of the Hadamard one, within the 
ADM-Hamiltonian formalism, with the result = or equivalently A = — (below we 
shall keep the value of A unspecified). 

The present paper's lineage is the harmon.c-coordinates approach HHEiEili. Its 
goal is the completion of the 3PN dynamics of compact binaries — equations of motion, 
Lagrangian, conserved integrals — in the frame where the center-of-mass is located at the 
origin of the coordinates. Our motivation is that the center-of-mass equations of motion 
constitute the needed starting point in many applications like the one in Ref. [2^. In 
Section II we recall the expression derived in Ref. |2^] for the position of the center-of-mass 
in an arbitrary harmonic-coordinates frame. In Section HI the individual positions of the 
particles in the center-of-mass frame are obtained as functions of the relative separation and 
velocity. We then compute the 3PN-accurate center-of-mass equations of motion. These 
equations are substantially simpler than in a general frame — though they are still quite 
lengthy (that is unavoidable at such a high post-Newtonian order). In particular we recover 
the center-of-mass equations of motion at the 2.5PN order derived by Lincoln and Will on 
the basis of the general-frame 2.5PN equations of Damour and Deruelle 0,0, IB]- In Section 
IV the 3PN relative Lagrangian (in harmonic coordinates), describing the conservative part 
of the dynamics, is obtained. Note that the center-of-mass relative Lagrangian does not 
strai ghtf orwardly follow from the general-frame Lagrangian of de Andrade, Blanchet and 
Faye (20] , because one is not a priori allowed to use in a Lagrangian some expressions which 
are consequences of the equations of motion derived from that Lagrangian. We found it 
convenient to derive the center-of-mass Lagrangian ab initio using some guess-work (i.e. 
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adjusting a set of coefficients in order to reproduce after Lagrangian variation the correct 
equations of motion). From that center-of-mass Lagrangian we then obtain in a standard way 
the Noetherian conserved energy and angular momentum, thereby completing our harmonic- 
coordinates approach. 

Further investigations are proposed. Section V deals with the connection between the 
center-of-mass Lagrangian and the center-of-mass ADM-Hamiltonian. We check that the 
center-of-mass reduction of the contact transformation worked out in Ref. between the 
harmonic-coordinates Lagrangian and the ADM-coordinates Hamiltonian is identical — as 
it must surely be — to the contact transformation connecting the center-of-mass versions 
of these Lagrangian and Hamiltonian. In the process we recover the 3PN Hamiltonian for 



the relative motion as computed by Jaranowski and Schafer Finally in Section VI 

we consider the problem of the stability against linear perturbations of circular orbits. We 
undertake the problem by perturbing both the equations of motion in harmonic coordinates, 
and the Hamiltonian equations in ADM coordinates (the two methods give equivalent re- 
sults). We obtain a gauge-invariant criterion for the stability of circular orbits up to the 
3PN order. 



II. THE CENTER-OF-MASS VECTOR POSITION 

Our study starts with the expression, derived in Ref. [2^, for the position G* of the 
binary's center of mass. In this Section we briefly review the construction of G*. Note that 
the center-of-mass position can also be interpreted as the gravitational mass-type dipole 
moment. Actually, using a slight abuse of language, by center-of-mass position G* we really 
mean the gravitational dipole (its dimension is that of a mass times a length). In a future 
work [2^ we shall show that the gravitational mass-type dipole moment which follows from 
a 3PN wave-generation formalism (instead of being inferred from the 3PN equations of 
motion) is in complete agreement with the present center-of-mass vector G*. 

By equations of binary motion (in a given coordinate system) we mean the acceleration 
a\{t) = dv\/dt of body A, where A = 1,2 and the spatial index i = 1,2,3, as a function 
of the positions and coordinate velocities f^(t) = dy'^^/dt. Eq. (7.16) in Ref. [l3| gives 
the 3PN equations of motion in the harmonic coordinate system. The 3PN Lagrangian in 
harmonic coordinates (considering only the conservative part of the dynamics) is given by 
Eq. (4.1) in Ref. Q; it takes the form 



L = L^[y^,\A] + ^^ipnIya, va] + ^^2PN[yA, va, a^] + — LspwiyA, va, a^] . (2.1) 

The successive post-Newtonian orders depend on the positions and velocities, and also, 
starting from the 2PN order, on the accelerations. The fact that a harmonic-coordinates 
Lagrangian necessarily becomes a "generalized" one (depending on accelerations) at the 
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2PN a ppr oximation has been proved by Damour and Deruelle j^^. At the 3PN order we 
found j20| that the Lagrangian also depends on accelerations, but it is notable that these 
accelerations are sufficient (i.e. there is no need to include derivatives of accelerations). 
Furthermore the dependence upon the accelerations at both the 2PN and 3PN orders is 
linear. Indeed one can always eliminate from a generalized Lagrangian, taking the form of 
a perturbative post-Newtonian expansion, a non-linear — for instance quadratic — term 
in the accelerations by adding a "double- zero" counter term, whose Lagrangian variation is 
zero on-shell and therefore which does not contribute to the dynamics (we refer to ^ for a 
general discussion on acceleration-dependent terms in a post-Newtonian Lagrangian). The 
conservative part of the equations of motion of body A (neglecting the 2.5PN and 3.5PN 
radiation damping terms) can be written with the help of the variational derivative of the 
Lagrangian as 

SL ^ dL _ d f dL\ d^ f dL 

It is important to remember that here the equations of motion are supposed to have been 
"order-reduced" using the equations themselves at lower post-Newtonian order (i.e. any a\ 
occuring in a post-Newtonian term must be replaced by its expression in terms of the ?/^'s 
and f^'s following the equations of motion). 

The existence of a center-of-mass integral for the 3PN dynamics is the consequence 
of its invariance under Lorentz transformations or boosts. The Lorentz invariance of the 
(harmonic-coordinates) equations of motion was established in Ref. jl^- Technically this 
means that a specific variant, defined in of the Hadamard regularization that respects 
the Lorentz invariance, is to be implemented. Consider an infinitesimal deformation of the 
path of the two particles, say Sy\(t) = y'\{t) —y\{t). Then the corresponding perturbation 
of the Lagrangian, i.e. 5L = L[-y'^, v^, a^] — L[yA, ^a, sla], reads, to the linearized order. 

It involves the total time derivative of the function 

Q = E (PA^yA + <1a^^a) , (2.4) 

A 

which is defined in terms of the momenta conjugate to the velocities and accelerations, 




5v\ dv\ dt yda'j^ 



^ This is consistent with a general argument of Martin and Sanz |2fll | that in coordinate systems which 
preserve the Lorentz invariance (as the harmonic coordinates do) the equations of motion at 2PN and 
higher orders cannot be derived from an ordinary Lagrangian. 
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Eqs. ()2.3p - ()2.5|) are nothing but the Noetherian equations in the case of a generahzed 
Lagrangian. 

In the case of a Lorentz transformation, the change in the position of particle A, to hnear 
order in the (constant) boost velocity V^, is 

Sy\ = -VH + ^V'y\v\ + ^(V^) . (2.6) 

Because the 3PN dynamics is invariant under Lorentz boosts, the change in the Lagrangian 
given by Eq. ()2.3|) must take the form of a total time derivative on-shell, i.e. when the 
equations of motion ()2.2|1 are satisfied. Hence there should exist a certain functional of 
the positions, velocities and accelerations such that (on-shell) 

SL = + 0{V^) . (2.7) 

Using this, together with the particular form of the transformation law ()2.(jj) . into Eq. ()2.3|1 
we readily obtain the conservation (on-shell) of the Noetherian integral = G^ — tP\ where 
P*, the total linear momentum, and G*, the center-of-mass position, are given by 



P^ = J2p\, (2.8a) 

A 

= + E {-^A + ^ Wa< + yWA< + v\q\v^^] ) . (2.8b) 

Since P* is itself constant [indeed, apply Eqs. ()2.3|) - ()2.4|) to the case of a constant spatial 
translation : 5y\ = e*], 

we find that the conservation of implies that 

f^p-. (2.10) 

[We neglect terms of order (9(c~®).] The center-of-mass vector G* is conserved in the rest 
frame where P* = 0; it will be zero, by definition, in the center-of-mass frame. 

Applying these considerations to the 3PN equations of motion and Lagrangian in har- 
monic coordinates, we found [2^ that indeed the variation of the Lagrangian uniquely defines 
some function (this is a confirmation of the boost symmetry of the equations of motion). 
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and from it we explicitly determined the center-of-mass vector position in an arbitrary 
harmonic-coordinates frame^ : 



mi 



+ 



+ 



mi 



+- 



mi 



m,2 



m2 

7 7 
- ^(^12-^1) - ^(^12-^2) ]v\ 



5 mim2 7m| ?)v\ 
4 rl^ Arl^ 8 

+ - ^{ni2Vlf - ^(^12^1) (77,12^2) + ^(^12^^2)^ 



235 mim2 235 m\ 

24 ri2 ' 24 ri2 

(53 3 
— {ni2Vif + -(7112^1)^(7112^2) + -{ni2Vi){ni2V2f 

5 15 . 1 

+— (?^i2^^2)^ - —{ni2Vi)vl - {ni2V2)vl + -(rzi2t'i)(t'iW2) 
iz o 4 



1 15 

+-:{ni2V2){viV2) - {ni2Vi)vl - —{ni2V2)v. 
4 o 



+ 



5vf 



+ —( ^{nuViY + J (7112^1)^ (7112^2) + :;^(7Zi2Vi)^(7Zi2V2)^ 

ri2 \ io o ib 

1 15 

+7(ni2t'i)(ni2t'2)^ - — (ni2f2)^ - — (ni2t'i)^fi 
4 lb lb 

X 11 53 3 

--{'>T'i2Vi){ni2V2)vl - —ini2V2yvf + —v^ + -(7112^1)^(^1^2) 
z o ib o 

3 5 

+ ^{ni2Vi){ni2V2){ViV2) + -{ni2V2Y{ViV2) 

17 1 5 

-5'yi('yi'y2) + -^-(^1^2)^ - ^(ni2'yi)^'y| - -(7112^1) (7112^2)^1 

+Y^{ni2V2) V2 + Y^ViV2 - y (^^i^^2)^^2 - Y^V2 



All-over this paper the gravitational constant is set to G = 1. 
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^?2 



7 29 2 
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3 
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~^ rf2 
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^ ~ V 1260 
_ 14351 
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40 , , 139 
3 V 1 2; 24 
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In 



r 



12 



-I ± 



r'r, J 



+ 1 ^ 2 + C 



(2.11) 



To the terms given above we must add those corresponding to the relabelhng 1^2. We 
denote by ri2 = |yi — y2|, = {y\ — y^^/ri2 and v\2 = v\ — v\ t^e relative particles' 
separation, unit direction and velocity. 

The expression (12. lip has been systematically order-reduced using the equations of motion 
and therefore depends only on the positions and velocities (no accelerations). Notice the 
appearance at the 3PN order of some logarithmic terms, containing two constants and 
(one for each body) having the dimension of a length. It was proved in Refs. 0, 



20| that these logarithms, and the r^'s therein, can be removed by an infinitesimal gauge 
transformation at the 3PN order. Thus we can refer to the r^'s as some "gauge constants", 
since they are merely associated with a choice of coordinate system, and thereby do not carry 
any physical meaning : they will always cancel out when deriving some physical, gauge- 
invariant, results. On the other hand we notice that is free of the physical regularization 
ambiguity A present in the equations of motion and Lagrangian. 

The previous derivation of the center of mass neglected the effect of radiation reaction. 
To take into account this effect we introduce some appropriate modifications at the 2.5PN 
order of the linear momentum and center of mass position : 



G' 



Am1m2 

5c^r^2 
4mim2 
5c5 



-n 



12 



2mi 



-"12 



+ 1^2 



Vf2 - 



ri2 

2(mi + 1712) 
ri2 



(2.12a) 
(2.12b) 



With these definitions, we find that the conservation laws ()2.9p and ()2.10|) . but now when 
taking into account the radiation-reaction effect, take in fact exactly the same form : 
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O - , (2.13a) 



dt 

f^P-.o(i). (2.13b, 

This finding is quite normal : recall that the total linear momentum of an isolated system 
is conserved up to the 3PN order included. Indeed the integral over the system of the local 
radiation reaction forces is a total time derivative at the 2.5PN order, and therefore it does 
not contribute to any change in the total linear momentum. The modification of the linear 
momentum by radiation reaction, or net radiation "recoil" of the source, is a smaller effect, 
of order 3.5PN — negligible in Eqs. ()2.13j) . 

III. EQUATIONS OF MOTION IN THE CENTER-OF-MASS FRAME 

The positions and velocities of the two particles in the center-of-mass frame at the 3PN 
order are obtained by solving the equation 

&[yA,VA] = oQ^^ , (3.1) 

where is defined in the previous Section. Obviously the solution must be determined 
iteratively, in a post-Newtonian perturbative sense, with systematic order-reduction of the 
equations. At the Newtonian order we get 



y\ = X^x' + ^ (^) ' (3.2a) 

yi = -X,x' + ^ (^) • (^-^^^ 
In this paper we employ the following notation. The relative binary's separation is 



x' = y\-yi, (3.3a) 

r = Ixl and = — . (3.3b) 

r 

For the relative velocity and acceleration we pose 

dx^ 

= — — = v\ — V2 and r = n ■ V , (3.4a) 

(Jjb 



(3.4b) 
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[r and were formerly denoted ri2 and v\2 in Eq. ()2.1H) ]. Concerning the mass parameters 
we denote 



mi 

— and 

m 

mi + m2 , 
mim2 



-^2 — , 

m 



m^ 



X1X2 and /i 



mi/ 



(3.5a) 
(3.5b) 
(3.5c) 



All the expressions that are written in the center-of-mass frame are conveniently param- 
eterized by m and the very useful mass ratio z/ (such that u = ^ in the equal-mass case 
and u in the test-mass limit for one of the bodies). Often it is convenient to consider 
reduced quantities, i.e. quantities divided by the reduced mass /x. 

The Newtonian solution ()3.2|1 is inserted into the IPN terms of Eq. (13.111 and we then 
obtain an equation for the IPN corrections in yl and 1/2 • Solving that equation we plug 
the result back into the IPN and 2PN terms of Eq. ()3.1|) and obtain the 2PN corrections 
in the same way. The process continues at the next order and this finally results in the 
3PN-accurate relationship between the individual center-of-mass positions y\ and 1/2 and 
the relative position and velocity v\ In the course of the computation we use for the 
order-reduction the center-of-mass equations of relative motion at the 2PN order — that is, 
at one post-Newtonian order before the 3PN order we want to reach. Since we give below 
the result for the 3PN equations, we do not detail this step and simply present the final 
expressions. They are in the form 



Vi 

2/2 



X2 + z/(Xi - X2)V + z/(Xi - X2)Qv' + O 



1 



Xi + z/(Xi - X2)V + z/(Xi -X2)Qv' + 



(3.6a) 
(3.6b) 



where all the post-Newtonian corrections, beyond the Newtonian result ()3.2|) . are propor- 
tional to the mass ratio v and the mass difference X1—X2. The two dimensionless coefficients 
V and Q depend on the mass parameters m, z/, the distance r, the relative velocity v'^ = 
and the radial velocity r = n ■ v : 



V 



V m 
3^;^ 3z/f^ 



m 



P 3r^iy 19v^ 3iyv^ 



r \ » 4 
2 u 



^m ^ 
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S = -r 



1 



16 



+ - 



m r 



16 



16 

2 ^ 16 



5r^f^ 21f^i/v^ 
+ 



16 

^7 4 
( UV 



16 

15z/2t;4 



+ - 



m 



7mr 



7P TSr^z/ 



3 

14351 



+ 4r^ z/^ + 



101 33z/i;2 



12 



1260 ^ 8 ~ T 



5 5r 
'5^2 I9r2j^ 



15^2 21i/w^ 



12 



24 

235 21 u 



+ 



24 



+ 3 z/2 



(3.7a) 



(3.7b) 



Jp to the 2.5PN order we find agreement with the circular-orbit limit of Eqs. (6.4) in Ref. 

30| (notice that the 2.5PN radiation-reaction term itself is proportional to the velocity and 
so it enters only the coefficient Q). 

In Eq. ()3.7|1 we find that the logarithms appear at the 3PN order and only in the 
coefficient V. They contain a particular combination Tq of the two gauge-constants r[ and 
r2 that is defined by 



(Xi - Xs) In r'^ = Xl In r[ - X| In . 



(3.8) 



This constant Tq happens to be different from a similar constant Tq which will have to be 
introduced to the 3PN equations of relative motion and Lagrangian [see Eq. 1)3.111) below]. 

The 3PN center-of-mass equations of motion are obtained in a straightforward way by 
replacing in the general- frame 3PN equations derived in Ref. [l^ (see Eq. (7.16) there) 
the positions by Eqs. ()3.6|) - ()3.7|) . and the velocities by the derivatives of Eqs. ()3.6p - ()3.7|) 
(applying as usual the order-reduction of all accelerations where necessary). Actually for this 
purpose we do not need the Eqs. ()3.6j) - ()3.7j) with the full 3PN precision; the 2PN- accurate 
ones are sufficient. We write the relative acceleration in the center-of-mass frame in the form 



dt r 

and find that the coefficients A and B are 



m 



1 + A) n' + Bv' 



+ 



(3.9) 
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A ^ — + v^ + 3uv^ 4 + 2i/ 

{ 2 r 

1 (15f^u 45r^i/2 9^2^^2 
\ 8 8 2 

+ — (-2 - 25 i/ - 2 i/^ - l^i^ + 2u^v'^ 
r \ 2 

/ 87 z/ 
+ ^ 9 + 



+ 



1 f 24 rz/?;^m 136 fz/m^ 



5 r 15 



1 r 35r6i/ 175f6i/2 i75f6i/3 15f^i/^;2 
^ c6\ 16~^ 16 16 ^ 2 

135 ^2 255 f^i/3^;2 15 r^^^;^ 237 f^^^^^ 
4 ^ 8 2 ^ 8 

— 

+ f 79 ^4 ^ _ 69rS^ - 30 - 121 ^ + 16 
r \ 2 

+ 20r2i/^'u2H hSu^v^ -lOu^v^ 

4 

/ 2 32573 llf2j,2 g;^5^2j,^2 26987 i/^;^ 
H r H ^ 7r u -\ 

3 2 _ USuTT^v^ -n0r'uln(^)+22iyvHn(^ 
64 Vr^/ Vr(, 

/ 41911Z/ 44Ai/ 71 i/^ Hun'^ 



+ — -16 + + } , (3.10a) 

V 420 3 2 16 ' ' ^ ' 



B = ^1 -4r + 2ri/| 



+ ^i^^ + 3r^z/2 ^ 2ri/%2 

c* 1^ 2 2 

m f 41 ru 
+ — 2f + + Afv 

r \ 2 

1 f 8 f 2 m 24 z/ 



1^ 5 r 5 
+ — < h 15 H h 12 r-^ i/ f ^ 



lllr'^z/2f2 657^1/^"^ .^.94 ^ . 4 

4 8 

m ( 329 V 59 r'"^ z/^ , a „ o 

H — \ h 18r'^z/ - 157^1/^2 - 27r z/^v^ - lOr z/ 

r V 6 2 



.3^2 



12 



/ 18169 fz/ . o 123rz/7r2 

+ — -4 r h 25 r + 8 r z/^ 

r2 V 840 32 



+44fz/ In (^^ H I . (3.10b) 



Up to the 2.5PN order the resuh agrees with the one given by Lincoln and Will j27|. At 
the 3PN order we have some gauge-dependent logarithms containing a constant Tq — dis- 
tinct from Tg introduced in Eq. ()3.8|) ^ — which is the "logarithmic barycenter" of the two 
constants r[ and : 

lnr[, = Xilnr'i + X2lnr^ . (3.11) 

In addition there is the physical ambiguity A due to the Hadamard self-field regularization 
(A cannot be removed by any coordinate transformation); it appears at the 3PN order in 
the ^-coefficient. 



IV. LAGRANGIAN AND NOETHERIAN CONSERVED INTEGRALS 

The Lagrangian for the relative center-of-mass motion is obtained from the 3PN center- 
of-mass equations of motion ()3.9j) - ()3.10j) in which one ignores for a moment the radiation- 
reaction 2.5PN term. The Lagrangian in harmonic coordinates will necessarily be a gen- 
eralized one, depending on accelerations, from the 2PN order. At the 3PN order, further 
acceleration terms are necessary but we do not need to include derivatives of accelerations. 
Furthermore we can always restrict ourselves to a Lagrangian that is linear in the accel- 
erations. Hence, our center-of-mass Lagrangian, denoted with a calligraphic letter C, to 
distinguish it from the general- frame Lagrangian L, is of the form 

£ = £n[x, v] + -^/:ipn[x, v] + -^/:2Pn[x, V, a] + -^£3pn[x, v, a] . (4.1) 

We recall that there is a large freedom for choosing a Lagrangian because we can always 
add to it the total time derivative of an arbitrary function. As a matter of convenience, we 
shall choose below a particular center-of-mass Lagrangian in such a way that it is "close" 
(in the sense that many coefficients are identical) to some "fictitious" Lagrangian that is 
obtained from the general-frame one given in Ref. 0] by the mere Newtonian replacements 
y\ — > y\ — s> —Xix\ We immediately point out that such a fictitious Lagrangian is not 

the correct Lagrangian for describing the center-of-mass relative motion. Indeed, the actual 
center-of-mass variables involve many post-Newtonian corrections given by Eqs. (j3.6|) - (j3.7|) . 



^ They are related by 
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so the actual center-of-mass Lagrangian must contain some extra terms in addition to those 

of the previous fictitious Lagrangian. However, we find that these extra terms arise only 
at the 2PN order and not before. We did not try to find a general method for obtaining 
systematically the center-of-mass Lagrangian given the general-frame one. Though such a 
method might exist it was in fact simpler to proceed by guess-work, i.e. to introduce some 
unknown coefficients in front of all possible types of terms, and to adjust these coefficients so 
that the Lagrangian reproduces the correct equations of motion. Our result (when divided 
by the reduced mass n = mi') is then 
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Witness the acceleration terms present at the 2PN and 3PN orders : our notation is a„ = a-n 
and = a-\ for the scalar products between a* = dv^/dt and the direction and velocity 
f The conservative part of the equations of motion is then identical (after order-reduction 
of the accelerations) to 

5x' ~ dxi dt \dv' J ^ dt^ \da' 7 ' 
From the Lagrangian one deduces the conserved energy and angular momentum using 
the generalized formulas [neglecting 0{c~^)] 



E = vY + a'q' - C , (4.4a) 
r = ei.k^xY' + v'q") , (4.4b) 

(the first one being a generalized version of the Legendre transform), where the conjugate 
momenta read 

6C ^dC d fdC 

« = &5 " • <^-^^' 

Alternatively one can compute the center-of-mass energy and angular momentum directly 
from the general-frame quantities E and J* given by Eqs. (4.2) and (4.4) in Ref. Q by 
replacing all variables by their center-of-mass expressions given by Eqs. ()3.6|l - ()3.7j] : the 
result is the same. For the energy we obtain 
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As for the center-of-mass angular momentum we get 
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(4.7) 



(The energy involves the regularization-ambiguity A, while the angular momentum is free 
of any physical ambiguity.) These quantities are conserved in the sense that their time 
variation equals the radiation- reaction effect. One can therefore modify them with terms of 
"odd order" to take into account the radiation reaction due to gravitational wave emission. 
For instance, in the leading 2.5PN radiation reaction one conventionally chooses that the 
right-hand side of the balance equations for energy and angular momentum take the standard 
form appropriate to the quadrupolar approximation. We then pose, 
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(4.8a) 
(4.8b) 



This choice is in agreement with the results of j^, 22, 23 1. Then we can easily check that 



dr 2 ■■■ (I 



(4.9a) 
(4.9b) 



where the Newtonian trace-free quadrupole moment is Qij = fi{x''x^ — ^6^^r^ 
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V. LAGRANGIAN AND HAMILTONIAN IN ADM-COORDINATES 



In Ref. ji^ (see also Ref. 0]) we determined the "contact" transformation between the 
particles' variables in harmonic coordinates and those in ADM (or rather ADM-type^) coor- 
dinates. By contact transformation we mean the relation between the particles' trajectories 
in harmonic coordinates, y\{t), and the corresponding ones in ADM coordinates, say Y\{t). 
We recall that in the contact transformation, i.e. 

Sy\{t) = YXit) - y\{t) , (5.1) 

the time coordinate t is to be viewed as a "dummy" variable^. There is a unique transforma- 
tion (15.11) such that the 3PN harmonic-coordinates Lagrangian of de Andrade, Blanchet and 
Faye j20] (in a general frame) is changed into another Lagrangian whose Legendre transform 
coincides with the 3PN ADM-coordinates Hamiltonian derived by Damour, Jaranowski and 
Schafer The explicit expression of this general-frame contact transformation can be 

found in Section 4.2 of Ref. j20|. 

Now we are in a position to obtain the relation between the relative separation vector 
= y\— y\ in harmonic coordinates and the one X* = Y"/ — Y2 in ADM coordinates (do not 
confuse the relative distances and X* between the two particles with the spatial position 
vector of some field event in these coordinates). Namely, 

5x' = X' -x' = 5y\ - 6yi , (5.2) 

where 5yl and 5y2 are given by Eqs. (4.8)-(4.10) in Ref. js^]. [We shall always view such 
equalities as ()5.2|1 as functional equalities, i.e. valid for any dummy time variable t.] One 
replaces in Eq. (15. 2|) all the variables by their center-of-mass counterparts following Eqs. 
fl3.6|) - ()3.7|) . Actually, since the contact transformation is already of relative order 2PN, the 
calculation is quite immediate and requires only the equations ()3.(jj) - ()3.7|l to IPN order. As 
a result the spatial separation vectors X* and in both coordinates (each one living within 
the spatial slice of constant time appropriate to each of the coordinate systems) are related 
to each other by 



m 



' ■ 1 
r V 



5 z/f m 
8 ^ 7 



3z/ 



H —v^ 



^ Strictly speaking, the ADM coordinates we are considering differ from the actual ADM coordinates at the 

3PN order by a shift of phase-space coordinates that is given in Ref. 13]. 
® The contact transformation is not a coordinate transformation between the spatial vectors in both 

coordinates, but takes also into account the fact that the time coordinate changes as well : i.e. 
= CiVA) — ^°(j/a)i'a/c, where £,'^{yA) denotes the four-dimensional change between the harmonic 

and ADM coordinates, when evaluated at the position ua — {t,y\). 
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(5.3) 



Below we shall deduce from this formula the radius of the circular orbit in ADM coordinates, 
say i?o, versus the one in harmonic-coordinates, i.e. ro; see Eq. ()6.19p . 

Now we look for the center-of-mass Lagrangian in ADM coordinates. Since in ADM 
coordinates the Lagrangian is "ordinary" (no accelerations) the contact transformation must 
be such that it removes the acceleration terms present in harmonic coordinates — more 
precisely, it must make them in the form of a total time derivative which is irrelevant to 
the dynamics. Following an investigation similar to the one in Section 3.2 of Ref. ji^l [see 
notably Eq. (3.18) there] we know that differs from C by two terms : (1) the functional 

variation of C induced at the linearized order by the contact transformation of the relative 
path as given by Eqs. ()5. 2^ - 1)5. 3|) : (2) the total time derivative of a function of the relative 
position and velocity. We can limit our consideration to the linearized order because 5x'^ is 
at least of order 2PN, so the non-linear terms do not contribute before the 4PN order and 
are negligible here. Hence we necessarily have the following functional equality (by which 
we mean the equality between functions of the same dummy variables x, v, a) : 



ADM 



X, vl = £ix, V, al 



5C ^ . dT 
ox* dt 



(5.4) 



in which is explicitly given by Eq. ()5.3|) . and where the function T is for the moment 
unknown^. We insist that in the present calculation the contact transformation 5a;* is known 
so that the only freedom left is the choice of T . This contrasts with our earlier study in Ref. 
|2o| where both the contact transformation of the individual paths, by\ and by\^ and some 
arbitrary function, say F, had to be varied and determined. The reason of course is that 
once 6y\ and known from Ref. 20| we have no choice for which must be equal 

to the center-of-mass reduction of the difference by\ — by^ [see Eq. (15. 2|) ]. Thus, despite 
the smaller freedom that we presently have in the adjustment of parameters, the calculation 
must work with that and not with another one. 



^ Note that the minus sign in front of the second term in Eq. 1)5.4(1 differs from the one in Eq. (3.18) of 
Ref. . The reason is because we have corrected a sign inconsistency in Ref. ■ namely the equation 
(3.12) there, together with the adopted definition fa ^ = y'j^ — y\, is inconsistent with (3.13) and (3.18); 
but this does not change any of the results of Ref. |20| . 
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The function T = J^[x, v] is not difficult to determine in order to match perfectly the 
ADM Hamiltonian. Notice that after adding the total time derivative of that function, not 
only has one been able to remove all the accelerations, but also one has gauged away all the 
logarithms which were present in the harmonic-coordinates Lagrangian. We get 
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Next, Eq. (j5.4p together with the explicit expressions (j5.3p and (|5.5p give the ADM center- 
of-mass Lagrangian. Once the calculation is done we have to express it using the names 
appropriate to the ADM variables : X* = + 6x^, which means the separation distance R, 
the relative square velocity V"^, and the radial velocity i? = N ■ V. The formula is 
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This is an ordinary Lagrangian and we apply the ordinary Legendre transform to obtain the 
Hamiltonian, which is a function of the conjugate momentum 



p2 = p2 and Pr = N-P . (5.7b) 
We find perfect agreement with the center-of-mass Hamiltonian derived in Refs. (ill . 0, : 
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^ -lecall that A is related to the so-called "static" regularization-ambiguity constant uo^ of Refs. 
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We have = according to Ref. 



the "kinetic" ambiguity constant tUk of Refs. [ll|, |l^ has been fixed to the value 



0. On 



the other hand, 



41 



the explicit Lorentz invariance of the equations of motion in harmonic coordinates |l6f il?! 
and, equivalently, by the requirement of existence of generators for the Poincare algebra in 
the ADM-Hamiltonian formalism jl3| . 

Finally let us present, for completeness, the formulas for the center-of-mass positions 
which are analogous to Eqs. ()3.6|l - ()3.7j) but in ADM coordinates. We have 



Yi = X2 + v{Xi - X2)V X' + u{Xi -X2)QV' + 



Xi + z/(Xi - X2)V X' + u{Xi -X2)QV' + 0[- 



(5.9a) 
(5.9b) 



where the post-Newtonian coefficients V and Q are given by 
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(5.10a) 
(5.10b) 



At the 2PN order the result is identical with the one given by Wex in his Appendix A j^^]^- 
By differentiating Eqs. (j5.9|) with respect to time we obtain the center-of-mass velocities 
and V2 in terms of the relative position X* and velocity K*. We have checked that by 
replacing into the obtained relations the velocities VI and V2 by their expressions depending 
on the conjugate momenta PI and P2 as deduced from the variation of the general-frame 
Lagrangian, and by expressing the velocity in terms of P* following the variation of the 
center-of-mass Lagrangian [Eq. (jS.Taj) ]. with both replacements being made with the full 
3PN accuracy, one ends up with the simple equations 



Pi = P' 



-PI 



2 ' 



(5.11) 



which are indeed the ones appropriate to a linear momentum that is conserved. 



VI. ON THE DYNAMICAL STABILITY OF CIRCULAR ORBITS 

As an application of the previous formalism let us investigate the problem of the stability, 
against dynamical perturbations, of circular orbits at the 3PN order. We propose to use 
two different methods, one based on a perturbation at the level of the equations of motion 
()3.9|l - ()3.1()j) in harmonic coordinates, the other one consisting of perturbing the Hamiltonian 
equations in ADM coordinates for the Hamiltonian ()5.8|) . We shall find a criterion for the 
stability of orbits and shall present it in an invariant way (the same in different coordinate 
systems). We shall check that our two methods agree on the result. 

We deal first with the perturbation of the equations of motion, following Kidder, Will 
and Wiseman i^] (see their Section III. A). We introduce polar coordinates (r, yj) in the 
orbital plane and pose u = r and uj = ip (beware that in this paper u = r, and not another 
standard notation in central force problems, u = 1/r). Then Eq. ()3.9j) yields the system of 
equations 



r = u , (6.1a) 



Since the contact transformation we consider relates together the conservative parts of the dynamics in 
harmonic and ADM-type coordinates (and affects only the 2PN and 3PN orders), the radiation-reaction 
damping term at the 2.5PN order in Eq. I)5.f Ob|l is the same as in harmonic coordinates. This is merely 
a definition of a particular ADM-type dynamics (a priori different from the one in ADM coordinates 
stricto-sensu), in which the "odd" terms, associated with radiation reaction, are the same as in harmonic 
coordinates. 
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+ ruj^, (6.1b) 

(6.1c) 

where A and B are given by Eqs. ()3.10p as functions of r, u and u (through f ^ = + r^cj^). 
In the case of an orbit which is circular apart from the adiabatic inspiral at the 2.5PN order 
(we neglect the 2.5PN radiation-reaction effect), we have r = ii = uj = hence u = 0. Eq. 
()6.1b|) gives thereby the angular velocity uq of the circular orbit as 
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(1 + A) . 



(6.2) 



Solving iteratively this relation at the 3PN order using the equations of motion ()3.9p - ()3.1U|) 
we obtain uq as a function of the circular-orbit radius tq in harmonic coordinates (the result 
agrees with the one of Refs. 0, E3|) • 
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[Please do not confuse the circular-orbit radius tq with the constant Tq entering the logarithm 

rTD.i 



at the 3PN order and which is defined by Eq. 

Now we investigate the equations of linear perturbations around the circular orbit defined 
by the constants ro, mq = [actually, if we were to include the radiation- reaction damping, 
uq = 0{c~^)] and ujq. We pose 
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ro + Sr , 



(6.4a) 
(6.4b) 
(6.4c) 



where 6r, 6u and Su denote some perturbations of the circular orbit. Then a system of 
linear equations follows : 



Sr = Su , 

6u = ao6r + [3q 5uj , 
5u = 7o 5m , 



(6.5a) 
(6.5b) 
(6.5c) 
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where the coefficients, which solely depend on the unperturbed circular orbit, read js^l 
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(6.6c) 



In obtaining Eqs. ()6.6p we use the fact that ^ is a function of the square through 
f ^ = + r'^uj'^, so that dA/du is proportional to u and thus vanishes in the unperturbed 
configuration (because u = 6u). On the other hand, since the radiation reaction is neglected, 
B also is proportional to u [see Eq. f)3.10b|) ]. so only dB/du can contribute at the zeroth 
perturbative order. Now by examining the fate of perturbations that are proportional to 
some e*'^*, we arrive at the condition for the frequency a of the perturbation to be real, and 
hence for stable circular orbits to exist, as being j^^l 



Co = -ao - /5o 7o > . 



(6.7) 



Substituting into this A and B at the 3PN order we then arrive at the orbital-stability 
criterion 
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(6.8) 



where we recall that tq is the radius of the orbit in harmonic coordinates. 

Our second method is to use the Hamiltonian equations based on the 3PN Hamiltonian 
in ADM coordinates given by Eq. ()5.8|) . We introduce the polar coordinates (-R, \1/) in the 
orbital plane — we assume that the orbital plane is equatorial, given by 6 = | in the 
spherical coordinate system (i?, 9, \1/) — and make the substitution 



= P, 



R 



p2 
i?2 



(6.9) 



into the Hamiltonian. This yields a "reduced" Hamiltonian that is a function of R, Pr and 
Pip : Ti = 7i \R, Pr, P.j,] , and describes the motion in polar coordinates in the orbital plane 
(henceforth we denote TC = TC^^^/ fi) . The Hamiltonian equations then read 
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(6.10a) 
(6.10b) 
(6.10c) 
(6.10d) 



Evidently the constant P,j, is nothing but the conserved angular-momentum integral. For 
circular orbits we have R = Rq {a constant) and Pr = 0, so 



f^[^o,0,P°]=0 



(6.11) 



which gives the angular momentum P^ of the circular orbit as a function of Rq, and 

(6.12) 

which yields the angular frequency of the circular orbit ujq — the same as in Eq. ()6.3p — in 
terms of -Rq ^ : 
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(6.13) 



We consider now a perturbation of the circular orbit defined by 



Pr = SPr , (6.14a) 

P^ = P^ + 6P^ , (6.14b) 

R = Ro + 5R, (6.14c) 

u = uo + Su . (6.14d) 



The last equation, 

oPr 

which is equivalent to R ~ const — Rq, is automatically verified because 7i is a quadratic function of Pr 
and hence dTL/dPR is zero for circular orbits. 
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It is easy to verify that the Hamiltonian equations ()6.1Up , when worked out at the hnearized 
order, read as 



6'R 

5uj 



— TTo 5R — Po 6P^i, 
0, 

o-Q 5Pr , 

Po 5R + To 6Pxs, , 



where the coefficients, which depend on the unperturbed orbit, are given by 



[-Ro, 0, P^] 



(^0 = -7^, [Ro, 0, P^] , 



To 



dPn' 

dp^^ 



[Ro, 0, P^] 



(6.15a) 
(6.15b) 
(6.15c) 
(6.15d) 



(6.16a) 
(6.16b) 
(6.16c) 
(6.16d) 



By looking to solutions proportional to some e*'^* one obtains some real frequencies, and 
therefore one finds stable circular orbits, if and only if 



Co = VToO-Q > 

Using the Hamiltonian ()5.8j) we readily obtain 



(6.17) 



Co 



ml m 

Ml ^ RW^ 



-9 + z/) + 



m 
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61 + 



135403 
1680 



4 

325 



-vr 



—A 
3 



(6.18) 



This result does not look the same as our previous result ()6.8|1 . but this is simply due to 
the fact that it depends on the ADM radial separation Ro instead of the harmonic one ro. 
Fortunately we have derived in Section V all the material needed to connect Ro to tq with 
the 3PN accuracy. Indeed, with Eqs. (|5.2|) - (|5.3p we have the relation valid for general orbits 
between the separation vectors in both coordinate systems. Specializing that relation to 
circular orbits we readily find 
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ro< 1 + 



+ O 




29 , 



m 



3163 21 2 

\ TT^ 

1680 32 
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In 



(6.19) 



The difference between Rq and tq is made out of 2PN and 3PN terms only. Inserting 
Eq. ()6.19|) into Eq. ()6.18|) and re-expanding to 3PN order we find that indeed our basic 
stabihty-criterion function Co comes out the same with our two methods. 

Finally let us give to the function C*o an invariant meaning by expressing it with the help of 
the orbital frequency c^o of the circular orbit, or, more conveniently, of the frequency-related 
parameter^" 



r3 j 



This allows us to write the criterion for stability as Co > 0, where Co 
gauge-invariant form (the same in all coordinate systems) 



(6.20) 
Cq admits the 



Co = 1 - 6xo + 14z/Xo + 



5954 123 



35 



16 



-vr 



44A 



14z/2 xl^O{xl) . (6.21) 



This form is more interesting than the coordinate-dependent expressions ()6.8|) or ()6.18|) . not 
only because of its invariant form, but also because as we see the IPN term yields exactly 
the Schwarzschild result that the innermost stable circular orbit or ISCO of a test particle 
(i.e. in the limit u ^ 0) is located at xisco = 1/6. Thus we find that, at the IPN order, 
but for any mass ratio u, 



X 



IPN 
ISCO 



1 

6 



(6.22) 



From the inverse of Eq. H6.3|l we obtain tq as a function of xq. For completeness we give the relations 
linking both and i?o to the XQ-parameter : 
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One could have expected that some deviations of the order of v aheady occur at the IPN 
order, but it turns out that only from the 2PN order does one find the occurence of some 
non-Schwarzschildian corrections proportional to v. At the 2PN order we obtain 



For equal masses this gives a;j™o — 0.187. Notice also that the effect of the finite mass 
corrections is to increase the frequency of the ISCO with respect to the Schwarzschild result 
(i.e. to make it more inward) : x^qq = | [l + ^z/ + C(z/^)]. Finally, at the 3PN order, for 
equal masses v = \ and for the value of the ambiguity parameter A = — (equivalent 
to tUs = 0), we find that according to our criterion all the circular orbits are stable. More 
generally, we find that at the 3PN order all orbits are stable when the mass ratio \^ v > 
where — 0.183. 

Note that the above stability criterion Cq gives an innermost stable circular orbit, when it 
exists, that is not necessarily the same as — and actually differs from — the innermost circu- 
lar orbit or ICO, which is defined by the point at which the center-of-mass binding energy of 
the binary for circular orbits reaches its minimum [33|. In this respect the present formalism, 
which is based on systematic post-Newtonian expansions (without using post-Newtonian re- 
summation techniques like Pade approximants jM differs from some "Schwarzschild-like" 
methods such as the effective-one-body approach 2^ in which the ICO happens to be also 
an innermost stable circular orbit or ISCO. 

As a final comment, let us note that the use of a truncated post-Newtonian series such 
as Eq. ()(j.2H) to determine the ISCO is a •priori meaningful only if we are able to bound 
the neglected error terms. Furthermore, since we are dealing with a stability criterion, it is 
not completely clear that the higher-order post-Newtonian correction terms, even if they are 
numerically small, will not change qualitatively the response of the orbit to the dynamical 
perturbation. This is indeed a problem, which cannot be answered rigorously with the 
present formalism. However, in the regime of the ISCO (when it exists), we have seen that 
is rather small: Xq — 0.2 (this is also approximately the value for the ICO computed in 
Ref. 01), which indicates that the neglected terms in the truncated series ()(j.2H) should not 
contribute very much, because they involve at least a factor Xg — 0.002. On the other hand, 
we pointed out that in the limit z/ — the criterion Cq gives back the correct exact result, 
\. This contrasts with the gauge-dependent power series ()6.8|) or ()6.18|) which give 



X 



ISCO 



only some approximate results. Based on these observations, we feel that it is reasonable 
to expect that the gauge-invariant stability criterion defined by Eq. ()fi.2Hl is physically 
meaningful. 
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